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ABSTRACT 



O 
O 

We have used the Keck I telescope to image at 11.7 /iin and 17.9 fim the 
dust emission around ( Lep, a main sequence A-type star at 21.5 pc from the 
Sun with an infrared excess. The excess is at most marginally resolved at 17.9 
/iin. The dust distance from the star is probably < 6 AU, although some dust 
may extend to 9 AU. The mass of observed dust is ~ 10 22 g. Since the lifetime 
of dust particles is about 10 4 yr because of the Poynting-Robertson effect, we 
robustly estimate at least 4 x 10 26 g must reside in parent bodies which may be 
(Sj , asteroids if the system is in a steady state and has an age of ~300 Myr. This 

mass is approximately 200 times that contained within the main asteroid belt 
in our solar system. 
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Subject headings: circumstellar matter — planetary systems — stars: individual 



1. INTRODUCTION 

Vega-like systems, first discovered by IRAS, are main sequence stars surrounded by 
dust. The operation of the Poynting-Robertson effect on micron-sized grains, responsible 
for much of the observed infrared excess, requires that the lifetime of these particles be 
significantly shorter than the age of the stars (Backman & Paresce 1993). Thus, the dust 
grains must be replenished from a reservoir such as collisions between larger bodies or the 
sublimation of comets. Studies of Vega-like systems thus afford the opportunity to study 
the evolution of large solids such as planets, comets, and asteroids. 

The identification of dust around main sequence A-type stars is usually made from the 
IRAS colors (e.g. Mannings & Barlow 1998; Silverstone 2000). Identifying objects with 
12 /im excess is difficult because the photosphere usually dominates the total flux at this 
wavelength. The bulk of the dust associated with main sequence A-type stars typically 
has temperatures T gr ~100 K, and semi-major axes >50 AU. For example, imaging of 
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HR 4796A (AOV, distance from Earth 67.1 pc) in thermal infrared (Koerner et al. 1998; 
Jayawardhana et al. 1998) and in near infrared scattered light (Schneider et al. 1999), 
revealed a dust ring with radius 70 AU, a distance scale similar to that of the Sun's Kuiper 
Belt. 

Dust grains located ~5 AU from a main sequence A-type star have T gr w 300 K, 
producing strong 12 /zm excesses. Searches for 12 /im excesses around main sequence stars 
have revealed that 300 K dust around A-K dwarfs is rare. In a survey of 548 A-K dwarfs, 
Aumann & Probst (1991) were able to identify IRAS 12 jum excesses only with f3 Pictoris 
and C Lep (= HR 1998 = HD 38678). Similarly, Cote (1987) derived unusually high 
temperatures for the dust around ( Lep compared to other Vega-like main sequence stars 
from IRAS 12, 25 and 60 /xm observations. Although coronagraphic images of (3 Pic (Kalas 
& Jewitt 1995) have revealed a massive dust disk which extends to a radius of almost 
1000 AU, no extended disk has been discovered around ( Lep. We have carried out a high 
resolution mid-infrared imaging study to learn about the dust properties, and the origin 
and evolution of the dust around ( Lep. 

C Lep is a main sequence A3Vn star with a Hipparcos distance of 21.5 pc and a 
fractional dust luminosity, L/^/L* = 1.7xl0~ 4 (Fajardo-Acosta, Telesco, & Knacke 1998). 
It has been identified as a member of the Castor moving group, which also contains 
Fomalhaut and Vega, based upon its kinematic properties (Barrado y Navascues 1998). 
Barrado y Navascues (1998) estimates its age to be between 100-400 Myr from the Lithium 
abundance of later-type stars in the Castor moving group and from the space motion of the 
group. Recently, Song et al. (2001) estimated an age of 50-350 Myr using Stromgren uvby/3 
photometry corrected for the effects of rapid rotation, to determine stellar parameters 
(effective temperature T e ff, surface gravity g and metallicity) which then are fit to stellar 
evolutionary tracks of Schaller et al. (1992). Lachaume et al. (1999) fit the position of ( 
Lep on the H-R digram with theoretical isochrones of Bertelli et al. (1994) and estimate a 
slightly older age of 200-500 Myr. 

2. OBSERVATIONS 

Our data were obtained on the night of 2001 February 5 (UT) at the Keck I telescope 
using the Long Wavelength Spectrometer (LWS), which was built by a team led by B. 
Jones and is described on the Keck website. The LWS uses a 128x128 SiAs BIB array with 
a scale at the Keck telescope of 0.08" pixel -1 and a total field of view of 10.2" x 10.2". We 
used the "chop-nod" mode of observing and two different filters: 11.2-12.2 /im and 16.9-18.9 
fim. We used Capella (a Aur) for flux and point-spread function calibrations. The data 
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were reduced at UCLA using standard LWS routines. 

We flux calibrate our data using the results for Capella that i 7 ^ (11.5 /im) = 176.1 
Jy and F^(19.3 /im) = 67.5 Jy (Gezari et al. 1987) extrapolated to our bands assuming 
that F v oc v 2 between 10 and 20 /im. For ( Lep, we find F^(11.7 /mi) = 1.60±0.20 Jy and 
Fj,(17.9 /im) = 0.93±0.07 Jy in agreement with previous observations by Fajardo-Acosta et 
al. (1998). We conservatively estimate the uncertainties associated with our measurements 
from the drift in the flux of Capella before and after our observations of ( Lep. Since less 
time elapsed between our observations of the standard star at 17.9 /im, the percentage 
errors are smaller at this wavelength. 

Determining the infrared excess of ( Lep requires subtracting a model for the stellar 
photosphere from the observed fluxes. We model the spectral energy distribution (SED) of 
( Lep with a 1993 Kurucz stellar atmosphere, assuming solar metallicity ([Z/H]=0.0). We 
assume no interstellar extinction because the star is only 21.5 pc from the Sun. We assume 
negligible circumstellar extinction since Li R jL if ~ 10~ 4 and we have no evidence in our 
data for an edge-on system. We find a best fit (minimum % 2 ) for the following parameters: 
stellar effective temperature T e // = 8500 K and surface gravity logg = 4.5 (see Figure 1). 
With this fit, we find photospheric fluxes of F 1/ (11.7 /im) = 1.24±0.01 Jy and F„(17.9 /im) 
= 0.53±0.01 Jy; thus, we find 11.7 /im and 17.9 /im excesses of 0.36±0.2 Jy and 0.40±0.07 
Jy respectively. 

We convolve models for the dust distribution with the point spread function (PSF) 
of Capella at 17.9 /xm to estimate the maximum distance of the dust from the star. 
Observations of Capella made 40 minutes before, 20 minutes after and 100 minutes after 
our observations of ( Lep varied somewhat. We use the PSF measured closest in time to 
our observations of ( Lep to construct our analysis. We assume the the dust is confined to 
a face-on ring, with an average distance of 0.10", 0.24" or 0.40" from the star and a width 
ranging between 0.20" and 0.24", and that 55% of the power is emitted by the point source 
and 45% of the power is emitted by the ring, consistent with our model for the photosphere 
and the total flux measured from ( Lep. We find a best fit for the model with the smallest 
dust disk (average radius 0.10"). The models including larger disks do not fit the surface 
brightness for the central pixel or the surface brightnesses for the outer pixels well. Thus, 
we estimate a maximum dust distance of 6 AU from the star. When models for the dust 
distribution are convolved with the PSF from earlier or later in the evening, the dust could 
extend as far as 9 AU from the star. 



3. DUST PROPERTIES 
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3.1. Minimum Grain Distance 

The minimum grain distance can be constrained from the temperature of the grains 
assuming that they are black bodies. We estimate a grain temperature of 320 K from the 
ratio of the 11.7 /iin excess to the 17.9 /im excess. Black bodies in radiative equilibrium 
with a stellar source are located a distance 

D= 1 -(^fR Sf (1) 

from the central star (Jura et al. 1998), where T* and are the effective temperature 
of the stellar photosphere and the stellar radius. We estimate the stellar luminosity from 
the bolometric magnitude using the Hipparcos V-band magnitude (my = 3.55 mag) and 
distance (21.5 pc) and a bolometric correction (Flower 1996) corresponding to an effective 
temperature T e ff = 8500 K. For ( Lep, we find a stellar luminosity = 14 Lq. The stellar 
radius is therefore /?* = 1.7 Rq. From equation (1) and the stellar properties summarized 
in Table 1, we find a minimum grain distance of 2.8 AU. 

We can additionally infer the dust temperature and distance by fitting a black body 
to the mid infrared fluxes reported for wavelengths longer than 10 /im. We find a best fit 
(minimum x 2 ) for the photosphere subtracted fluxes for T^ust — 230 K and D = 5.4 AU, 
consistent with our observations. 



3.2. Circumstellar Dust Grain Size 

A lower limit to the size of dust grains orbiting a star can be found by balancing the 
force due to radiation pressure with the force due to gravity. For small grains with radius 
a, the force due to radiation pressure overcomes gravity for: 

a < 3L*Q pr /(167rGM*cp s ) (2) 

(Artymowicz 1988) where and M* are the stellar luminosity and mass, Q pr is the 
radiation pressure coupling coefficient and p s is the density of an individual grain. Since 
radiation from an A- type star is dominated by optical and ultraviolet light, we expect that 
2na/\ ^> 1 and therefore the effective cross section of the grains can be approximated by 
their geometric cross section so Q pr ~ 1. Based upon T e ff and L*, the inferred stellar mass 
is 1.9 Mq (Siess, Dufour, & Forestini 2000). With p s = 2.5 g cm -3 , the minimum radius 
for grains orbiting ( Lep is a = 1.7 /im. 

We can estimate the average size of the grains assuming a size distribution for the 
dust grains. Analogous to the main asteroid belt in our solar system and as expected from 
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equilibrium between production and destruction of objects through collisions (Greenberg & 
Nolan 1989), we assume 

n(a)da = n Q a~ p da (3) 

with p ~ 3.5 (Binzel, Hanner, & Steel 2000). If we assume a minimum grain radius of 
1.7 /mi, we find an average grain radius < a > = 2.8 /im, if we weight by the number of 
particles. 



3.3. Mass of Circumstellar Dust Around £ Lep 

We can estimate the minimum mass of dust around ( Lep assuming that the particles 
have < a > ~ 2.8 /jm; if the grains are larger, then our estimate is a lower bound. If we 
assume a thin shell of dust at distance, D, from the star and if the particles are spheres of 
radius, a, and if the cross section of the particles equals their geometric cross section, then 
the mass of dust is 

M d > f^PsD 2 <a> (4) 

(Jura et al. 1995) where L IR is the luminosity of the dust. If D = 6 AU, p s = 2.5 g cm -3 , 
and < a > = 2.8 pm, then M d = 1.6 xlO 22 g. 



3.4. Lifetime of Circumstellar Grains 

One mechanism which may remove particles is Poynting-Robertson drag. The 
Poynting- Robertson lifetime of grains in a circular orbit, a distance D from a star is 

/4tt < a > p s \ c 2 D 2 
tpR = —r~ (5) 



(Burns et al. 1979). With the parameters given above and = 14 Lq, the Poynting- 
Robertson lifetime of the grains is tpR = 1.3 xlO 4 years. Since this timescale is significantly 
shorter than the stellar age (t age ), we hypothesize that the grains are replenished through 
collisions between larger bodies. By analogy with the solar system, we propose that the 
parent bodies are asteroids. 
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3.5. Mass of the Parent Bodies 

We can estimate the total mass contained in parent bodies around ( Lep assuming a 
steady state. If M PB denotes the mass in parent bodies, then we may write 

M PB > ^t age (6) 

With M d = 1.6 x 10 22 g, t PR = 1.3 x 10 4 yr and t age = 300 x 10 6 yr, we estimate a total 
mass of parent bodies ~ 4 x 10 26 g, approximately 200 times the mass of the main asteroid 
belt in our solar system (Binzel et al. 2000). From equations (4) and (5), we find 

Mp B > ^SSl ( 7 ) 

Since L IR is well measured, our estimate for the lower limit on the mass of parent bodies, 
Mpb, is better constrained than our estimate for the mass of observed dust. 



4. DISCUSSION 

C Lep is distinctive because the dust around it is warm and close to the star compared 
to the dust in other well known debris disk systems: (3 Pictoris, HR 4796A, Vega, and 
Fomalhaut. We can quantify this difference by comparing the ratios of 60 /iin excess to the 
25 /xm excess. For ( Lep, this ratio is at least a factor of 2 smaller than for any of the other 
stars considered (see Table 2). 

The lack of dust at distances greater than 6 AU raises the possibility of the presence 
of a planet, sculping the disk and confining the dust. Such a planet could increase the 
eccentricity of the orbits of the putative asteroids and thus drive them into mutual collisions 
to produce the observed dust. 



5. CONCLUSIONS 

We have obtained high resolution mid infrared images of ( Lep at 11.7 fim and 17.9 
jum using the LWS on the Keck I telescope. 

1. ( Lep possess a strong mid infrared excess which is at most marginally resolved 
and thus surrounding dust probably lies within 6 AU although some dust may extend as 
far as 9 AU. This result is consistent with the temperature of the grains derived from the 
mid-infrared photometry. 
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2. Since the estimated Poynting-Robertson lifetime of grains with < a > = 2.7 /im 
is 1.3 x 10 4 yr, significantly shorter than the age of ( Lep, the grains must be replenished 
from a reservoir such as collisions between larger bodies, perhaps asteroids. For ( Lep, we 
robustly estimate that the minimum mass of parent bodies is ~ 4 x 10 26 g assuming that 
the system is in a steady state. This mass is approximately 200 times the mass of the Sun's 
main asteroid belt. 

This work has been supported by funding from NASA. We thank M. Sykes, A. 
Weinberger, and B. Zuckerman for their comments. 
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Table 1. £ Lcp Properties 


Quantity 


Adopted Value 


Reference 


Spectral Type 


A3Vn 


1 


Distance 


21.5 pc 


2 


Effective Temperature (T e //) 


8500 K 




Surface Gravity (log g) 


4.5 




Stellar Radius (R#) 


1.7 RO 




Stellar Luminosity (L») 


14 L Q 




Rotational Velocity {vsmi) 


202 km/sec 
1.7xl0~ 4 


1 


Fractional Dust Luminosity 










Age 


50-500 Myr 


3, 4, 5 



References. — (1) HofHeit & Warren 1991; (2) Hipparcos; (3) Barrado y Navacues 1998; (4) Lachaume et 
al. 1999; (5) Song et al. 2001 



Table 2. Disk Properties 



Star 


Photosphere+Excess a 

F v (60^™) 
F l ,(25fj,m) 


Excess b 

F v (60^m) 
F v (25tJ,m) 


T color 

(K) 


C Lep 


0.32 


0.30 


370 d 


/3 Pic 


2.2 


2.3 


100 


HR 4796A 


2.2 


1.8 


110 


Vega 


0.85 


4.4 


80 


Fomalhaut 


2.0 


11 


70 



TRAS data without color corrections 

b IRAS data with color corrections and photosphcric subtraction 

c derived from F^(60 (im)/F u (25 /im) for the excess radiation given in the previous column 
d & black body plus stellar photosphere fit to all of the fluxes reported for wavelengths greater than 10 ^m 
yields a lower temperature, T gr = 230 K 
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-10 12 3 

LOG Wavelength (//m) 

Fig. 1. — Spectral energy distribution (SED) for £ Lep with ultra-violet fluxes from TD 1 
(Thompson et al. 1978) plotted as squares and IUE plotted as asterisks and UBV photometry 
(Cousins 1984) plotted as crosses. Previous JHK and mid-infrared photometry from the 
IRTF (Fajardo-Acosta et al. 1998) are shown with dark blue error bars, from the IRAF 
Faint Source Catalog (Moshir et al. 1989) with green error bars and from ISO (Habing et 
al. 2001) with red error bars. Our 11.7 /im and 17.9 fim photometry, as reported here, is 
shown with magenta error bars. Overlaid is a 1993 Kurucz model for a stellar atmosphere 
with T e ff = 8500 K, surface gravity log g = 4.5 and solar metallicity ([Z/H]=0.0) plotted to 
minimize x 2 . The dotted black line is a best fit to the mid infrared fluxes at wavelengths 
longer than 10 /im assuming a single temperature black body (T = 230 K) plus the stellar 
photosphere. 
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Fig. 2a.— Radial profile for C Lep (a(J2000.0)=05 /l 46 m 57 s , <5(2000.0) = -19 o 49'19") at 
11.7 /Ltm. The solid black line is the point spread function for the standard star Capella 
(a(J2000.0)=05H6 m 41 s , <5(2000.0)=+45°59'53") normalized to the observations of C Lep, 
while the error bars show the azimuthally averaged data for ( Lep. 
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Fig. 2b. — Radial profile for £ Lep at 17.9 //m. The solid black line is the point spread 
function for the standard star Capella, normalized to the observations of ( Lep, while the 
error bars show the azimuthally averaged data for ( Lep. The colored lines are best fit 
models (minimum x 2 ) f° r the spatial extent of the dust convolved with the point spread 
function of Capella. The red, green and blue models assume a central point source and a 
ring with average angular radius 0.10", 0.24" and 0.40" respectively. The \ 2 values for the 
different curves here are 4.5, 5.2, and 9.3. In each model, 55% of the power is emitted by 
the point source and 45% of the power is emitted by the ring, consistent with our model for 
the photosphere and the total flux emitted from ( Lep. 



